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We develop a theoretical description of electro-magnon solitons in a coupled ferroelectric- 
ferromagnetic heterostructure. The solitons are considered in the weakly nonlinear limit as a mod¬ 
ulation of plane waves corresponding to two, electric- and magnetic-like branches in the spectrum. 
Emphasis is put on magnetic-like envelope solitons that can be created by an alternating electric 
field. It is shown also that the magnetic pulses can be amplified by an electric field with a frequency 
close to the band edge of the magnetic branch. 
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Multiferroic materials, i.e., materials exhibiting cou¬ 
pled order parameters, are in the focus of current re¬ 
search. These systems offer not only new opportunities 
for applications but also provide a test ground for ad¬ 
dressing fundamental issues regarding the interplay be¬ 
tween electronic correlations, symmetry, and the interre¬ 
lation between magnetism and ferroelectricity []]-[2 . Here 
we address magnetoelectrics which possess a simultane¬ 
ous ferroelectric-magnetic response. A interesting aspect 
is the non-linear nature of the magnetoelectric excita¬ 
tion dynamics, which hints at the potential of these sys¬ 
tems for exploring nonlinear wave-localization phenom¬ 
ena, such as multicomponent solitons Bi , nonlinear 
band-gap transmission [6, 0, and the interplay between 
the nonlinearity and Anderson localization [8[. In this pa¬ 
per we aim at exciting robust magnetic signals by means 
of electric fields. Particularly, we consider a multifer¬ 
roic nano-heterostructure consisting of a ferromagnetic 
(FM) part deposited onto a ferroelectric (FE) substrate. 
As demonstrated experimentally, under favorable condi¬ 
tions, a coupling between the ferroelectric and the ferro¬ 
magnetic order parameters may emerge (this coupling is 
referred to as the magnetoelectric coupling), thus allow¬ 
ing one to control magnetism (ferroelectricity) by means 
of electric (magnetic) fields. Here we consider the case 
when the multiferroic structure is driven by an electric 
field with a frequency located within the band-gap of the 
FE branch and in the band of the magnetic-excitation 
branch. For a proper choice of the electric-field frequency 
(that follows from the electro-magnon soliton theory de¬ 
veloped below) it is possible to excite propagating mag¬ 
netic solitons. In addition, we point out a possibility for 
the amplification of weak magnetic signals, which sug¬ 
gests the design of a digital magnetic transistor, where 
the role of the pump is played by the electric field. 

Examples of the two-phase multiferroics under study 
[M2, are BaTiOs/CoFe 204 or PbZri_ x Ti x 03 /ferrites. 
The developed model will be applied to a system where 
the FE and FM regions are coupled at an interface whith 


a weak magnetoelectric coupling. The theory is, how¬ 
ever, more general and can, in principle, be applied to 
single-phase magnetoelectrics |l3l~ll5j| . For the creation 
of electro-magnon solitons, which is the subject of the 
present work, a two-phase multiferroic structure is more 
appropriate, as it allows to generate and manipulate iso¬ 
lated FE or FM signals away from the interface. 

Both single- and two-phase multiferroics may be mod¬ 
eled by a ladder consisting of two weakly coupled chains: 
One chain is ferroelectric (FE) built out of unidimen¬ 
sional electric dipole moments, P n . The second chain 
is ferromagnetic (FM), composed of classical three- 
dimensional magnetic moments, S n , where n numbers 
the site in the lattice. Each chain is characterized by an 
intrinsic nearest-neighbor coupling, and each P n is cou¬ 
pled to S n via inter-chain weak magnetoelectric coupling. 
For a discussion of the microscopic nature of this coupling 
we refer to [16] . We assume the direction of FE dipoles at 
some arbitrary angle with respect to FM anisotropy axis 
£, as depicted in Fig. la. The magnetoelectric coupling 
will cause a rearrangement of magnetic moments. Let the 
new ground-state ordering direction of FM be the axis z, 
and cj) is the angle between z and anisotropy axis £. The 
magnetic field h(t) is applied along z, and 6 is an angle 
between z and FE moments (see Fig. la). So (Po) stands 
for the FM (FE) equilibrium configuration. We will con¬ 
sider perturbations around the equilibrium. Defining the 
scaled dipolar deviations p n = (P n — Po)/Po and the 
scaled magnetic variables s n = S n /So, the Hamiltonian 
is written as 
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FIG. 1. (Color online) (a) A schematic of the multiferroic 
ladder built of FE and FM chains. Arrows indicate the di¬ 
rections and the magnitudes of the electric dipole moments 
and magnetic moments in the course of the soliton advance¬ 
ment along the ladder, (b) Mutual orientations of FE and 
FM ground-state vectors, in the framework of Hamiltonian 
©• (c) Dispersion relations for FE (upper) and FM (lower) 
branches of the multiferroic composite. The arrow indicates 
the selected carrier frequency. 


where H sp stands for the linearized interfacial magne- 
toelectric coupling between the FM and the FE chain 
[l6| . Hp is FE part of the energy functional for N- 
interacting FE dipole moments 171. fl8|. Further, Ao is 
a kinetic coefficient; aj is the nearest-neighbor coupling 
constant; a and f3 are second- and forth-order expansion 
coefficients of the Ginzburg-Landau-Devonshire (GLD) 
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potential SzL 0 near the equilibrium state Pq. H s 
stands the ferromagnetic contribution (20[, where J is 
the nearest-neighbor exchange coupling in the FM part. 
D i = DSo cos 2 (p and D 2 = DSo are anisotropy con¬ 
stants, and D is the uniaxial anisotropy constant along 
axis £ (see Fig. lb). 

We operate with dimensionless quantities by using the 
scaling t -A uuot with ujq = y/aj/a 0 (for the examples 
shown below ujo ~ 10 12 rad/sec). The other parameters 
of the model <7, A, /?, J, Z>i, D 2 are scaled with ujq. 
The scaled quantities are indicated by omitting the tilde 
superscript. The time evolution is governed by 

r)^ x 

= - J [ s n ( s n—1 + s n+l) “ s n ( s n-l + s n+l)] “ 

-2 DM, 

c) 

« J [4 (4- 1 + 4+1) - 4 (4-1 + 4+l)] + 

+ 2D lSn s n ~ gPnSn, ( 2 ) 

d 2 p n 

= ~ a P n ~ PPn + (Pn -1 - 2 Pn + Pn+ 1) + 9S 

As we are interested in small perturbations, s®, s ^ and 
p n are much less than unity and the approximate equality 
4 = 1 - (4) 2 / 2 - (4) 2 / 2 is justified. 

We seek weakly nonlinear harmonic solutions to Eq. 
©, with a frequency uj and a wavenumber /c, in the form 
of a column vector (s®, s^,p n ) = Rexp [i (uot — fcn)]+c.c., 
where R is a set of complex amplitudes R = (a, 6, c), 
and c.c. stands for the complex conjugate. Neglecting 
higher harmonics in Eq. ©, we find the set of nonlinear 
algebraic equations 


W*R=<2 n \ (3) 

where the matrix and source are, respectively, 


W = 


iu Jsin 2 (/c/2) + 2D 2 0 

[J sin 2 (/c/2) + 2 D 2 \ iuj g 

g 0 uj 2 — a — sin 2 (fc/2) 


/ b(\b\ 2 + |a| 2 )( J — J cos k + D 2 ) — b*(b 2 + a 2 )(Jcos k — J cos 2/c — \ 

_ j _ a (|5|2 ^ | a |2^ j _ j cog 7. ^ X ) + a*(6 2 + a 2 )(J cos — Jcos 2/c — Di) J. (4) 
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The linear limit amounts to the set of linear homoge¬ 
neous algebraic equations W * R = 0. The solvability 
condition Det ^ W^j = 0 leads to two branches of the dis¬ 
persion relation uj{k) which are shown in Fig. 1(c), with 
the corresponding amplitude set, R = (a, 6, c), where b 
and c are expressed via the arbitrary constant a: b = 
—iau;/ [ J sin 2 (/c/2) + D 2 ], c — gaj [a + sin 2 (/c/2) — u \. 


We call a dispersion branch ferroelectric (defining its fre¬ 
quency uje and labeling the amplitude with index E ), if 
it has \c\ > \a\ [the red curve in Fig. 1(c)], while a ferro¬ 
magnetic branch (ujm) is defined by the relation \c\ < \a\ 
(the blue curve in the same figure). 

Of a particular interest is the case when the system 
is excited at an edge (at the left one, for the sake of 
definiteness), with a frequency uj s which falls into the 
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bandgap of FE mode and, simultaneously, the propa¬ 
gation band of the FM one, as shown by the arrow 
in Fig. 1(c). The dispersion relation with the fixed 
frequency, uj = uj s becomes then a cubic equation for 
sin 2 (fc/2). For uj s belonging to the band of FM mode and 
bandgap of FE one, the cubic equation yields two com¬ 
plex wavenumbers, associated with FE and FM modes, 
and a real one, corresponding to the FM mode. These 
three wavenumbers determine a set of three orthogonal 
eigenvectors, R E = ( a E ,b E ,c E ), R M = 
and Rm = ( a M> c m); where the first two corre- 
spond to complex FE and FM wavenumbers, kE and 
k^, respectively, while the last one is related to the real 
wavenumber, k^- In linear systems, the solutions with 
the complex wavenumbers are evanescent waves localized 
at the left edge of the multiferroic chain. Thus, the solu¬ 
tion for the vector function, F n = (s®,s^,p n ), is 

= A{t) R E e <h, ‘ t - |fcB l n + c.c. 

F™“ = B(t)R^e iuJst -^ n + c.c. (5) 


where the amplitudes A(t) and B(t) may vary slowly in 
time. 

As mentioned above, the solutions corresponding to 
the complex wavenumbers are localized at the boundary. 
To examine the possibility of a solitonic self-localization 
of the third solution with a real wavenumber (cf. Refs, 
for similar solutions in multiferroic models), we 
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consider the nonlinear frequency shift produced by the 
small terms Q nl in Q. Assuming a shifted frequency, 
uj + Scj instead of uj, the matrix W is substituted by a 
modified one, W + 6 W, with the diagonal matrix JW = 
iScu • Diag(l, 1, 2u). We also define a row vector L = 
(a', b\ c') which solves for the equation L* W = 0. Then, 
multiplying both sides of Eq. m on L, we obtain the 
nonlinear plane-wave frequency shift: 


Slj (fc, |a| 2 ) 


(l, «•■)/(: 


L* JW*R 


( 6 ) 


From these results, operating with the en¬ 
velope function <p(n, t ) defined from F^ + = 
p(n, £)e z ( Wst_ ^ n ) (a^, 6^, c]^), one can derive 
the nonlinear Schrodinger equation (NLS), cf. Refs, 
mm in the form: 


2 i 


dp dp 

~di +V &n 


+ W "dn? + _ ^ 


which gives rise to the respective envelope-soliton solu¬ 
tion with the FM-like localized mode being written as 


F“* 


cosh 



vt ) 


a M \ 

b M I + c - c " (8) 

C M J 


The velocity, dispersion, and the nonlinearity coefficients 
are 

_ didM ,, _ d 2 (JM * _ d [< 5cdM(k , |a| 2 )] 

dk ’ ^ dk 2 ’ ^[| a l 2 ] 

(9) 

The carrier frequency of this soliton is defined by the 
dispersion relation [the lower blue curve in Fig. 1(c)]: 

U s = u M (k) + SuiM (k, |a| 2 ) /2. (10) 

Thus, one can generate both the FE evanescent (|5]) 
and FM solitonic ^ modes, driving the left edge of the 
chain at the same frequency, uj s . It is possible to produce 
a combination of these solutions too. Generally in nonlin¬ 
ear systems, linear combinations of particular solutions 
is not another solution but if the solutions are far sepa¬ 
rated, which makes interactions between them negligible, 
the linear combination 

Fn = /iF® + h F™~ + F^ + (11) 


is still a solution of the nonlinear problem. In the weakly 
nonlinear limit it is even possible to construct a solu¬ 
tion for the case when particular modes overlap (i.e., the 
magnetic soliton is located near the edge), adding a time- 
dependent phase to each term © and (|5]) in the sum [26[ . 
For instance, one can consider an approximate solution 
at the left edge of the ladder, n = 0, in the form of 


m ~ 


(t) 


F M+ e » M W 


F 0 = /iFjfe^O + / 2 F£ 

( 12 ) 

where, in the weakly-nonlinear limit, the phases T are 
proportional to the wave amplitudes. Hence the waves 
do not gain significant phase shifts due to interaction 
effects, if their relative group velocity is not negligible. 
In this case, all phase shifts may be neglected. 

Our particular aim is to create an FM soliton by ex¬ 
citing only the FE degree of freedom at the edge, i.e., 
Sq = Sq = 0- To this end, we choose A(t) = B(t) = 


sech 


a^vt 


A/8 w" 


seeking to impose the following 


vector relation at the edge, n = 0: 


(0,0, Po) = 


(Ji Re + / 2 R; 


’M 


■R-m) 


0 iu) s t 


cosh 


a M V 


ty/ A/8( 




T c.c. (13) 


Using now the orthogonality of eigenvectors R, we 
readily get the appropriate expression for po : 


Po(t) = 


D + 


( c m) cosh a^vty/ A/8 cj / 


T + C.C. 


(14) 


Further, it is possible to compute the coefficients f\ and 
/ 2 , using the same orthogonality property: 


| r m| m , |Hj^| ( c m) 

|Re| 2 (c+)*’ /2 |R M | 2 (c+)*' 


(15) 
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In numerical simulations, if we make po a function of 
time as per Eq. m, keeping the magnetic moments 
pinned at the boundary, it is possible to excite the FE 
and FM evanescent waves (|5j), and also propagating FM 
soliton dHJ. These simulations correspond to an experi¬ 
mental setup with pinned boundary conditions at both 
FE and FM edges, and to the application of the elec¬ 
tric field E(t) = po(t) according to Eq. (fl4l) at the first 
cell of the FE chain. In this way, one can realize the 
excitation of magnetic solitons in the FM chain of the 
multiferroic ladder via an electric (rather than magnetic) 
field by virtue of the magnetoelectric coupling. 

For an assessment of the above, we performed full nu¬ 
merical simulations with the following values of the nor¬ 
malized parameters 

a = 0.2; p = 0.1; J = 1; D x = 0.1; D 2 = 0.2; g = 0.1. 

__ (16) 

These values correspond to BaTiO^/Fe J27|, [28|. Fur¬ 
thermore, we assume for the FE second and fourth or¬ 
der potential coefficients au/(a EE ) = 2.77 • 10 7 [Vm/C], 
ch/(a EE ) = 1-7 • 10 8 [Vm 5 /C 3 ], and for the FE coupling 
coefficient dj/(a EE ) = 1.3 • 10 8 [Vm/C], the equilibrium 
polarization Po = 0.265 [C/m 2 ], and the coarse-grained 
FE cell size ap e = 1 [nm]. The FM exchange interac¬ 
tion strength is J = 3.15 • 10“ 20 [J], the FM anisotropy 
constant is D = 6.75- 10“ 21 [J], and the ME coupling 
strength is go ~ 10 -21 [Vm 2 ]. 

We drive the left boundary of FE according to m 
with a driving frequency uj s = 0.4 and an amplitude 
a M = 0-07 and apply pinned boundary conditions for 
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FIG. 2. (Color online) Numerical simulations of the creation 
of a ferromagnetic soliton by the electric field, using Eq. © 
with parameters (H6l) . (a,b) Show a comparison between the 
analytical (solid lines) and the numerical (points) results for 
the soliton’s spatial profile at different moments of time. The 
analytical solution is taken according to Eq. CD with the 
coefficients (H5l) . (c,d): The propagation of the FM soliton in 
the ferroelectric and the ferromagnetic layers, respectively. 



FIG. 3. (Color online) Results of numerical simulations of the 
amplification of a magnetic signal in the multiferroic chain (a). 
The space-time evolution of excitations in the ferroelectric 
part, (b) The amplified propagation of the magnetic soliton 
in the ferromagnetic part. The arrow indicates the moment 
of the injection of the magnetic signal. Inset (c) presents the 
time dependence of the energy of the input magnetic signal 
at the FM edge, n — 0, while the energy of output signal is 
measured at site n = 200. 


FM, .sg = = 0. The results are displayed in Fig. 2, 

where the comparison of numerical simulations and ap¬ 
proximate analytical solution (HD are shown at different 
moments of time [Figs. 2(a,b)]. 

Next, let us envisage the possibility of amplifying mag¬ 
netic pulses: We apply a continuous electric signal with 
the frequency uj s which is slightly below the FM band 
boundary, cjm( 0), keeping FM moments pinned at the 
edge. In this setting, and for small driving amplitudes, 
no energy is transmitted through the chain. Both FM 
and FE modes are evanescent and described by the solu¬ 
tions ©• A propagating FM soliton emerges only if the 
electric-field amplitude attains the band gap transmis¬ 
sion threshold [6|, [29|: This happens if the amplitude is 
large enough so that a solution of the nonlinear dispersion 
relation m for real wave number k exits. Then, if one 
keeps the amplitude of the electric field just slightly be¬ 
low this band gap threshold, a small-amplitude in-phase 
magnetic signal coupled to FM chain allows to pass the 
threshold. This gives rise to a large-amplitude FM soli¬ 
ton propagating through the multiferroic, see Fig. 3. In 
this way, one can realize an amplification of the magnetic 
pulses by electric field. It may happen that almost all the 
energy of the electric field will be transferred to the ferro- 
magnet chain, and the corresponding amplification rate 
may achieve values as much as ^ 100. For instance, in the 
simulations presented in Fig. 3 we choose the driving fre¬ 
quency and the amplitude of the electric field uj s = 0.229 
(that is « 30 GHz in real units) and (po) max = 0.414, 
while the FM signal amplitude is ( s o) max = 0-02. 
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In this paper we do not address dissipation effects 
which, in principle, could be taken into account by intro¬ 
ducing the conventional Landau-Lifshitz-Gilbert damp¬ 
ing term in the magnetic part of the evolution equations 
©, as well as damping terms in the electric part. Here, 
we assume that dissipation has no qualitative effects for 
the considered length and on the time scales compara¬ 
ble with magnetic/electric signal transmission (that is 
10 -9 sec) and do not consider thus the respective terms 
in the evolution equations. 

Concluding, an electro-magnon soliton theory is de¬ 
veloped and the results are applied for electric field- 
induced magnetic soliton generation. A proper choice 
of pump electric field parameters enables an amplifica¬ 
tion of magnetic signals. In the amplifying regime the 
total (pump+signal) amplitude overcomes the band-gap 
transmission threshold and the energy of electric field 
is completely transferred to the magnetic soliton. As we 
have shown above, substantial (more than 100 times) am¬ 
plification of the magnetic input/output signals could be 
realized. 
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